Dilatometric Study of LiHoF4 In a Transverse Magnetic Field 
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Theoretical and experimental work have not provided a consistent picture of the phase diagram of 
the nearly ideal Ising ferromagnet LiHoF4 in a transverse magnetic field. Using a newly fabricated 
capacitive dilatometer, we have investigated the thermal expansion and magnetostriction of LiHoF4 
in magnetic fields applied perpendicular to the Ising direction. Critical points for the ferromagnetic 
phase transition have been determined from both methods in the classical paramagnetic to ferro- 
magnetic regime. Excellent agreement has been found with existing experimental data suggesting 
that, in this regime, the current theoretical calculations have not entirely captured the physics of 
this interesting model system. 



I. INTRODUCTION 

The study of ideal Ising systems has attracted consid- 
erable attention, owing to the computational tract ability 
of theoretical models and the discovery of some materi- 
als which exhibit nearly ideal Ising characteristics^. One 
such material is LiHoF4, which becomes ferromagnetic 
at the Curie temperature, Tc = 1.53 k!^. The mag- 
netic behavior is due the spin on the Ho^+ ion. The 
Ising anisotropy is provided by the crystalline electric 
field which splits the multiplet of the Ho^+ states yield- 
ing a groundstate doublet with the spins pointing par- 
allel or antiparallel to the [001] axis of the crystal. The 
first excited state is separated from the groundstate by 
approximately 11 K^. As the temperature is lowered 
the magnetic ordering proceeds via a second-order phase 
transition. In principle, this is a well-understood clas- 
sical phenomenon driven by thermal fluctuation^. In 
this case, the system dimension is equal to the critical 
dimension leading to measurable logarithmic corrections 
to mean- field critical exponents^. 

The application of a magnetic field transverse to the 
Ising (easy) axis of LiHoF4 splits the ground state dou- 
blet, and introduces quantum mechanical tunneling be- 
tween the Ising states of the Ho^+ ions in the presence 
of the crystalline electric field. The associated quantum 
fluctuations lead to a suppression of the Curie temper- 
ature with increasing field^, culminating in a quantum 
phase transition when the transition temperature is sup- 
pressed to zero. This physics is thought to be captured by 
one of the simplest theoretical models to describe a quan- 
tum phase transition, namely the transverse-field Ising 
model. Originally proposed by de Gennes^, the Hamilto- 
nian for this model is given h^ 
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Here Jij are coupling constants between spins, a^'^'^ 
are Pauli operators, and F is the effective transverse field 



FIG. 1: A schematic phase diagram of the transverse field 
Ising model where T is the temperature and F is the effec- 
tive transverse field parameter. The solid line separates the 
magnetically ordered phase from the paramagnetic phase. 



parameter that represents the mixing of the two Ising 
states in the presence of the physical transverse magnetic 
field. For F 7^ 0, quantum mechanical fluctuations in the 
Ising order occur even at zero temperature. Increasing 
the fluctuations by increasing the transverse field even- 
tually leads to a disordering of the magnetically ordered 
phase at a critical value of F = Fc. This is the quan- 
tum critical point. At finite temperature, the quantum 
fluctuations lead to a suppression of the transition tem- 
perature compared with the zero-transverse-field value. 
A phase diagram is shown schematically in Fig. [l] 

Such quantum phase transitions have attracted consid- 
erable attention in the physics community recently^. In 
this respect, the transverse field Ising model provides one 
of the simplest realisations of such fascinating and topical 
physics. Paradoxically, for the application of the trans- 
verse field Ising model to LiHoF4, whilst the quantum 
critical regime appears moderately well-understood^ ^-^, 



a significant discrepancy exists between the experimental 
data and the theoretical description of the phase diagram 
in the predominantly classi cal reg ime close to the zero- 
field transition temperature^ ^^ ^^. It is in this regime 
where a complete description requires a full understand- 
ing of the effects of temperature on the quantum fluctua- 
tions. Ultimately, this may prove to be the most relevant 
physics for understanding the novel behaviour related to 
quantum fluctuations in real systems occurring at finite 
temperatures^^. 

Using measurements of magnetic susceptibility, Bitko 
et al.^ were the first to report the phase line of the ferro- 
magnetic state in LiHoF4 in a transverse magnetic field 
using measurements of magnetic susceptibility. A rea- 
sonable theoretical fit to the data was presented using 
mean- field theory albeit with two adjustable parameters; 
the Lande g-factor and effective dipole coupling strength 
to characterise the effect of magnetic field and tempera- 
ture, respectively. 

Recognizing the limitations of such a simplified theo- 
retical approach, Chakraborty et al}^ developed a full 
microscopic Hamiltonian that purportedly captured all 
the relevant physics of this transverse field Ising system, 
notably the effects of quantum fluctuations and the do- 
main structure of the ferromagnetic state. Parameters 
in the theory were fixed where possible by results from 
spectroscopic^^ and susceptibilit}!^ measurements leav- 
ing one adjustable parameter as the antiferromagnetic 
Heisenberg exchange coupling. This quantity was intro- 
duced to explain the measured zero-field value of the 
ferromagnetic ordering temperature which is lower by 
some 25% over the quantum Monte Carlo value fixed by 
the known parameters'^. Using quantum Monte Carlo 
numerical techniques, a phase line was computed that 
agreed well with the data particularly in the low temper- 
ature limit where quantum fluctuations are significant. 
Nonetheless, a discrepancy with the data close to the 
zero-field transition was apparent and reported as re- 
flecting the potential uncertainties in the crystal- field 
parameters for LiHoF4. 

Further theoretical work on the transverse field Ising 
model as it relates to LiHoF4 followed by R0nnow et 
al^. This study addressed neutron scattering spectra'^ 
as well as the phase line measured earliei^. The critical 
points determined from neutron scattering experiments 
lacked close agreement at low transverse field compared 
to previous susceptibility work^. For the theoretical cal- 
culation, the crystal-field p aramete rs were derived from 
spectroscopic measurements^^^Sn^ on Ho-doped LiYF4 
which agreed closely with previously determined values^. 
As in the Chakraborty worlff^, the computed phase line 
overestimated the rate of suppression of the ferromag- 
netic state at very low transverse fields. However, given 
that now different crystal-field parameters have led to a 
similar discrepancy between theory and experiment as 
had been seen earlier, it is unlikely that this is the origin. 

Most recently, a substantial theoretical effort by Tabei 
et al^ has been reported that, in particular, attempts 



to resolve the differences between the theoretical phase 
line and experiment at low transverse fields^^^. Us- 
ing the fact that in this regime the quantum fluctua- 
tions are small, they were introduced perturbatively into 
a classical Hamiltonian. The calculation then proceeded 
using classical Monte Carlo techniques that were both 
significantly different and simpler to implement than 
their quantum Monte Carlo counterparts. The possible 
sources of discrepancy were separated into ones which 
were computational in nature and ones that resulted from 
inadequacies in the model Hamiltonian. The authors 
concluded that the discrepancy was not of computational 
origin, nor did it stem from uncertainties in crystal-field 
parameters. Consequently, it was most likely due to 
shortcomings of the model Hamiltonian (Eq. IT]) used. 
Furthermore, they also made the suggestion that the ex- 
perimental determination of the phase-diagram should 
be revisited. This provided substantial motivation for 
the study reported here. 

As outlined above, there have been several attempts to 
improve the theoretical understanding of the phase dia- 
gram of LiHoF4, while little has been done to repeat the 
experimental measurement of the phase diagram. Here 
we report a detailed experimental study of the phase line 
using thermal expansion and magnetostriction measure- 
ments. We focus exclusively on the high temperature - 
low field regime where there are markedly few data points 
from the original experimental invest igatioiP and where 
the theory is most challenged. Our results show good 
agreement with the existing experimental data indicating 
that there is still more theoretical work to be done. The 
rest of this article will briefly discuss our experimental 
technique, the methods used to determine critical points 
herein, and then present the results from this investiga- 
tion into the transverse field phase diagram of LiHoF4. 



II. APPARATUS AND EXPERIMENTAL 
METHOD 

Previous experimental investigations have made use of 
magnetic susceptibilit}'^, neutron scattering'^ and spe- 
cific heaiP (zero-field only) to determine critical points 
in the phase diagram of LiHoF4. For the study reported 
here, capacitive dilatometry has been used to measure 
the thermal expansion and magnetostriction of a LiHoF4 
sample. A detailed description of this device will be lim- 
ited in this article, and we refer interested readers to a 
future publication'^. 

The LiHoF4 sample is taken from a commercially pro- 
duced single crystal'^. It is shaped as a roughly semi- 
circular plate, where the easy axis [001] of LiHoF4 is ori- 
ented approximately 25° from the flat edge of the semi- 
circle. The sample is approximately 1 mm thick, 5.5 
mm tall, and about 3.5 mm at it's widest point. The 
easy [001] axis of the LiHoF4 sample was determined to 
within ±1° using a commercial Laue diffraction setup, 
and the ends of the crystal subsequently polished to be 



orthogonal to this easy axis (within ±5°). When placed 
in the dilatometer, length changes along the [001] axis 
were measured as temperature and magnetic field were 
varied. 

The dilatometer used here is conceptually similar to 
that reported by Schmiedeshoff et al?'^ . It is fabricated 
from 99.995% pure silver and differs in that the moveable 
capacitor plate is mounted to a BeCu spring known for 
its fiexibility at low temperatures. The sample is then 
pressed against this spring, and does not make direct 
contact with the moveable capacitor plate. Samples are 
thermally linked to a sample mount (made of 99.995% 
pure silver) using a silver paint^^. An uncalibrated Cer- 
nox CX-1030 resistance thermometer is also fastened di- 
rectly to this mount and is used to verify the thermal 
stability of the sample. 

A pumped ^He cryostat is used to cool the sample 
mounted in the dilatometer to approximately 1 K. The 
temperature is controlled using a Lakeshore 331 Tem- 
perature controller and measured by a calibrated Cer- 
nox CX-1030 thermistor located on the cold stage of the 
pumped ^He cryostat. An experimental trial consists 
of a temperature sweep through the region of interest 
{^ 1.2 K to 1.8 K) at a rate of 7.5 mK/min. 

Capacitance was measured using a commercial 
Andeen-Hagerling 2500A capacitance bridge, with a res- 
olution of 10~^ pF, which equates to a dilation of 
0.002 A at 17 pF. The measured capacitance is con- 
verted to capacitor plate separation using an expression 
for a tilted plate capacitor^ 
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FIG. 2: Zero- field measurement of the thermal expansion co- 
efficient, a, of LiHoF4 and the cell effect as a function of 
temperature. The critical temperature, Tc is indicated, be- 
ing the first deviation from the paramagnetic state. Inset: 
Normalised thermal expansion coefficient and specific heat 
measurements^ as a function of reduced temperature. 
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where a^g is the thermal expansion coefficient of silve: 
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Where eo is the permittivity of free space, A is the ef- 
fective capacitor plate area, C^ax is the maximum mea- 
sured capacitance before the capacitor plates short, and 
C is the measured capacitance. Prior to beginning a 
measurement, the maximum capacitance (C^ax, typically 
^ 50 pF) is measured. This is the measured capacitance 
just before the plates touch and electrically short (leading 
to C = pF), and provides enough information to take 
into account the tilt of the capacitor plates. Once typ- 
ical operating conditions are taken into account, length 
changes greater than 0.2 A can be resolved. A measured 
change in the capacitor plate spacing is directly related 
to a change in the sample length, dD = —dL. However, 
the measured length changes will include the effect of the 
dilatometer as well as the LiHoF4 sample. To decouple 
the behavior of the sample from the measured capaci- 
tance changes two pieces of information are required: the 
behavior of the dilatometer with a silver sample (known 
as the cell effect), and the behavior of the dilatometer 
with LiHoF4. Ultimately, the quantity of interest is the 
thermal expansion coefficient, a(T)^ which is the nor- 
malised rate of change of length with temperature. The 
relation between the measured thermal expansion coef- 
ficient, and the thermal expansion coefficient of LiHoF4 



III. RESULTS & DISCUSSION 
A. Thermal Expansion 

Zero field: In Fig. [2J the thermal expansion of LiHoF4 
is plotted in the absence of any magnetic field. The tran- 
sition from the paramagnetic to the ferromagnetic state 
as the temperature is lowered results in a sudden jump in 
the thermal expansion coefficient. Consistent with pre- 
vious experimental^ and theoretical worlPl', the critical 
temperature, Tc^ is defined as the first deviation from the 
behaviour in the paramagnetic state as the temperature 
is decreased. This point is determined from the intersec- 
tion of a polynomial fit to the temperature dependence of 
the thermal expansion above the transition and a linear 
fit to the steepest section of the thermal expansion during 
the transition (increasing the order of either fit made no 
appreciable difference). The error in this determination 
was estimated by varying the range of the fitted curves 
near the critical region by up to ±0.05 K (10% of the 
temperature range), and observing the range of Tc's that 
result. 

Using the above analysis, the transition temperature 
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FIG. 3: Measurement of the thermal expansion coefficient, a, 
for LiHoF4 up to an applied transverse field of 1.25 T. The 
critical temperature at each field is indicated. The curves at 
different field strengths are offset by 3 • 10~^ per curve. 



FIG. 4: Measurement of the thermal expansion coefficient, a, 
for LiHoF4, from 1.25 T to 2.75 T. The critical temperature 
at each field is indicated. The curves at different fields are 
offset by —5 • 10~^ per curve. 



in zero field is T^ = 1.532 ± 0.005 K. In the inset to 
Fig. |2J we plot the thermal expansion data and high- 
resolution specific heat data^. One can relate the thermal 
expansion coefficient {a) to the specific heat {Cy) through 
the Griineisen parameter (7) and the bulk modulus {B): 
a(T) = (7/35) Cy(T). Assuming the Griineisen param- 
eter and the bulk modulus are temperature independent, 
it is apparent that the discontinuities observed in specific 
heat measurements should be reproduced in a thermal ex- 
pansion measurement. To aid close comparison, the ther- 
mal expansion data is inverted and both measurements 
are normalised to their peak values while being plotted 
against reduced temperature with respect to their indi- 
vidual transition temperatures. The qualitative agree- 
ment is excellent (as seen inset in Fig.^. Quantitatively, 
there is a smaller discrepancy between the absolute value 
of transition temperatures (1.532 ± 0.005 K and 1.5384 
K) while the transition width is larger by 4.4 mK in the 
thermal expansion measurements. Both effects can rea- 
sonably be attributed to the swept temperature method 
that was used in this experiment. 

Transverse Magnetic Field: The evolution of the tem- 
perature dependence of the thermal expansion in trans- 
verse field is shown in Figs. [3]and|4] In these figures the 
curves at each field have been offset for clarity. Three 
qualitative features are apparent in this data set. First, 
the temperature dependence of the thermal expansion 
coefficient in the paramagnetic state evolves from low 
positive gradient towards a negative gradient of slightly 
larger magnitude at a field of 1 T and then returns to a 
positive gradient by a field of 2 T. Second, in the ferro- 
magnetic phase, the temperature dependence transitions 
from concave downwards in zero field and low fields to 



concave upwards at fields above 1 T. Third, the transi- 
tion region itself is modified to accommodate these vari- 
ations, starting as a step downwards upon entering the 
ferromagnetic phase and evolving to being a step upwards 
for fields of 1 T and above. This is accompanied by a sig- 
nificant broadening of the transition region as the field is 
increased. 

As in the zero field measurement, the transition tem- 
perature is defined as the first deviation from the mono- 
tonic behaviour that is considered the paramagnetic 
phase. The same analysis as described above is used to 
extract a transition temperature at each magnetic field. 
As the transition width broadens with increased mag- 
netic field, the uncertainty in the transition temperature 
increases and assumes a larger error bar. 



B. Magnetostriction 

Magnetostriction, A = AL/L, is the change in length 
of a material in response to a magnetic field. The mag- 
netostriction of LiHoF4 was measured at T = 1.3, 1.35, 
1.4 and 1.8 K and the results are plotted in Fig. |5]as a 
function of B^ . Since the measurement at T = 1.8 K 
is entirely in the paramagnetic region, by subtracting 
this from the reponse at the lower temperatures, AA = 
A(T) — A(1.8K), we expose the magnetic response more 
clearly. The inset shows the first derivative of A A at 
1.3, 1.35 and 1.4 K with respect to the squared trans- 
verse field. While the raw magnetostriction signal does 
not exhibit any distinct features during the phase transi- 
tion as seen in the thermal expansion, there is a distinct 
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FIG. 5: Magnetostriction, A = ^, as a function of transverse 
field squared, B^) measurements performed on LiHoF4, in 
the minimum demagnetization orientation. Inset: ^^ as a 
function of B^ . The critical field is the determined when the 
derivative equals zero. 



FIG. 6: Measurement of a for LiHoF4 at ±2.5 T in the max- 
imum (solid lines) and minimum (dashed lines) demagnetiza- 
tion orientations. Crosses show the calculated critical tem- 
peratures. The red (blue) dotted line indicates the critical 
temperature determined from the maximum (minimum) de- 
magnetization orientation, identical within error. 



change in this derivative which we interpret as the field 
at which ferromagnetic state is entered as the field is low- 
ered. The critical field is therefore determined to be the 
point at which the first derivative equals zero. The un- 
certainty associated with this value is dominated by the 
uncertainty of the magnetoresistance of the thermistor 
used to control the temperature. 



C. Demagnetization Effects 

The irregular shape of the sample means that the de- 
magnetization factor due to uncompensated spins on the 
surface of the material is difficult to calculate. However, 
we can assess the impact of the demagnetization effect 
on our results by orienting either the semi-circular face 
or the thin edge of the crystal perpendicular to the trans- 
verse field. These two experimental orientations repre- 
sent a maximum and minimum demagnetization factor 
for this particular sample. In the following we estimate 
the difference in the magnitude of the demagnetization 
factor between the two orientations. 

The demagnetization correction can be written as 



n,. 



H^ 



Applied 



GM, where H is either the inter- 



nal or applied magnetic field, G is the geometric factor 
resulting from the crystal geometry, and M is the sample 
magnetization. To a first approximation, the difference 
between i^intemai and i^AppUed will be determined wholly 
by G, when the only experimental change is the orienta- 
tion of the crystal and we assume the magnetization to 
be the same in each situation. Beleggia et alM- proposed 
a scheme for magnetostatic mapping of shapes to equiva- 
lent ellipsoids whose geometric factors can be more easily 
calculated. Following this work, we can estimate the de- 



magnetization factors of the LiHoF4 sample used in this 
work by a rectangular prism with dimensions L^ = 1 mm, 
Ly — 2.1 mm, and L^ = 5.5 mm. 

Using the Mathematica code provided in the above ref- 
erence, we find the demagnetization factors G^ = G^ax = 
7.55, Gy = G^i, = 1.68, and G. = -0.32. Thus, the 
correction for the maximum demagnetization orientation 
will be approximately 4.5 times larger than in the min- 
imum demagnetization orientation. By repeating the 
same experiments using these two different sample orien- 
tations, we can compare demagnetization corrections. As 
expected, in all orientations, the observed zero-field tran- 
sitions are identical and occur at T^ = 1.532±0.005 K. 
Fig. [6] shows the results of such an investigation at ilf = 
2.5 T. As can be seen, a striking feature of this test is 
the qualitative difference in the temperature dependence 
of the thermal expansion between each orientation. The 
transition temperature which is extracted in each case is 
Tc(max) = 1.35±0.03 K and Tc(min) = 1.37±0.03 K for 
each orientation as shown by the crosses in the figure. 
Consequently, in spite of a difference in the demagneti- 
zation factor of almost five, there is only a small differ- 
ence in the transition temperature which is less than our 
uncertainty in each value. This suggests that the demag- 
netization correction will lead to changes in T^ that are 
at most within our current error estimates. 

The preceding analysis assumes the sample is single do- 
main. Calculation of the demagnetisation factor is a far 
more complicated issue in the presence of domainP^^^. 
Nonetheless, even in the presence of a magnetic domain 
structure, we would still anticpate an orientation de- 
pendence to the demagnetisation factor for this sample 
and the empirical evidence from this study is that this 



correction is small. Theoretical analysis of the domain 
structure in LiHoF4 has been carried out by Biltmo and 
Henelius^. We speculate that differences in the domain 
structure that could occur for each orientation may ex- 
plain the qualitative orientational difference in the signal 
that is measured. A better understanding of the micro- 
scopic origin of the thermal expansion coefficient and the 
behavior of domains in the region of the transition would 
be required to explain this observation. 



D. Torque Effects 

The potential misalignment of a ferromagnetic sample 
in a transverse magnetic field may lead to a torque on 
the sample which, if the sample were to rotate or bend, 
could lead to a signal indistinguishable from thermal ex- 
pansion or magnetostriction. In order to verify that any 
torque applied to LiHoF4 was not distorting our abil- 
ity to measure the phase transition, thermal expansion 
was measured in positive and negative transverse fields, 
while in both the minimum and maximum demagnetiza- 
tion orientations. 

As shown in Fig. [6J the measured thermal expansion 
for either field direction is identical. Consequently we 
conclude that there is not a significant rotation or bend- 
ing of the sample due to any putative torque. This is 
consistent with a small sample magnetic moment due to 
magnetic domains. 



IV. DISCUSSION & CONCLUSIONS 

Combining the ferromagnetic to paramagnetic criti- 
cal temperatures from thermal expansion measurements 
with critical fields from magnetostriction measurements 
produces the phase diagram shown in Fig. [7| The some- 
what broad error bars reflect the broadening of the tran- 
sition region and consequently the associated difficulty 
in assigning a transition temperature rather than a loss 
of precision in our technique. Also plotted are the criti- 
cal points from earlier magnetic susceptibilit}'^ and neu- 
tron scattering^ work, and from the latest numerical 
simulations^. There is excellent agreement between the 
experimental data from our technique and the magnetic 
susceptibility, thereby confirming the discrepancy with 
the current theoretical simulations. Qualitatively, what 
appears to occur is that when a transverse field is ap- 
plied to LiHoF4 in the ferromagnetic state the system is 
initially impervious to the field. It is only when the field 
rises above 1 T, which is approximately 25% of the field 
required to completely suppress the ferromagnetism at 
zero temperature, that the transition temperature starts 
to be suppressed. The origin of this imperviousness to 
the transverse field is unknown and would appear to be 
an aspect of the physics that is not captured by the trans- 
verse field Ising model (Eq. [l]) used to describe LiHoF4 
used in numerical simulationP^^IIIIl What is now im- 
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FIG. 7: Transverse field versus temperature phase diagram of 
LiHoF4, as measured by this study using thermal expansion 
(squares) and magnet ostr iction (circles) data. Also plotted is 
previous experimentaP^^ and theoretical Monte Carlo^ work. 
(Dotted line is a guide to the eye.) 



portant to assess is whether these results reveal a generic 
lack of understanding of the impact of quantum fluctua- 
tions at finite temperature, or the simpistic nature of the 
transverse field model as applied to LiHoF4. 

One aspect of physics that is specific to LiIIoF4 is the 
effect of magnetoelastic coupling on the antiferromag- 
netic or quadrupolar exchange parameters as discussed 
in earlier theoretical studies^ ^^. These couplings com- 
pete with the ferromagnetic coupling and can, in princi- 
ple, affect the transition temperature. R0nnow et al?^ 
estimated that a strain, 613 ~ 10~^ would be required to 
produce a noticable effect in the transition temperature 
due to the quadrupolar interactions. Our measurements 
of magnetostriction are only along the c-axis (633) so we 
are unable to address this issue directly. The overall mag- 
nitude of the magnetostriction measured in this study is 
Al/l ~ 10~^. Consequently, unless the magnitude of the 
in-plane magnetostriction is dramatically different (i.e 2 
orders of magnitude larger), it is unlikely that this could 
lead to a significant restoration of the critical tempera- 
ture to compensate for the reducing effect of the quantum 
fluctuations. Further studies of magnetostriction along 
the a-axis would be required to provide a more complete 
picture. In particular, a qualitative change in magne- 
tostrictive behavior between measurements shown here 
at temperatures very close to Tc(0) and those at lower 
temperatures would be illuminating. 

On a related point, we note that the magnetic field 
scale of 1 T also signifies a qualitative change in the be- 
haviour of the transition in thermal expansion between 
the paramagnetic and ferromagnetic state, as seen in 
Figs. |3|4[ In the absence of a detailed microscopic the- 
ory describing the features of the transition, we are not 
able to assess whether this is a related phenomena or 
merely coincidental. As already discussed, an experimen- 



tal investigation of domain wall motion to complement 
existing theoretical work^^ would be invaluable. Careful 
analysis of the evolution of the transition in any other 
technique may also provide additional insight. 

Finally, investigation of other Ising systems, for exam- 
ple Ho(OH)3 and Dy(0H)3, as proposed by Stasiak et al. 
(Ref . [28|) . could determine if this discrepancy between the 
theoretical model and experimental work is particular to 
this material, or is a generic feature of the effect of a 
transverse field on Ising magnetism in real systems. 

In conclusion, we have used measurements of thermal 
expansion and magnetostriction to map the phase dia- 
gram of LiHoF4 in a transverse magnetic field at tem- 
peratures close to the zero-field classical phase transition 
where quantum fiuctuations are small. The resultant 
phase line separating the paramagnetic phase from the 
ferromagnetic phase is in agreement with earlier exper- 
imental studies^ and disagrees with current theoretical 



results based on the transverse field Ising modeP^^. The 
crucial issue now is whether the disagreement pertains to 
details of the LiIIoF4 system or whether it is a generic 
feature signifying a lack of understanding of the effect of 
quantum fiuctuations at finite temperature. 
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